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ABSTRACT 

In this paper we present an error analysis for polarimetric data obtained with dual-beam instru- 
ments. After recalling the basic concepts, we introduce the analytical expressions for the uncertainties 
of polarization degree and angle. These are then compared with the results of Monte-Carlo simula- 
tions, which are also used to briefly discuss the statistical bias. Then we approach the problem of 
background subtraction and the errors introduced by a non-perfect Wollaston prism, flat-fielding and 
retarder plate defects. We finally investigate the effects of instrumental polarization and we propose a 
simple test to detect and characterize it. The application of this method to real VLT-FORS1 data has 
shown the presence of a spurious polarization, which is of the order of ~1.5% at the edges of the field 
of view. The cause of this effect has been identified with the presence of rather curved lenses in the 
collimator, combined with the non complete removal of reflections by the coatings. This problem is 
probably common to all focal-reducer instruments equipped with a polarimetric mode. An additional 
spurious and asymmetric polarization field, whose cause is still unclear, is visible in the B band. 
Subject headings: Instrumentation: polarimeters - Methods: data analysis 
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1. INTRODUCTION 

Performing polarimetry basically means measuring 
flux differences along different electric field oscillation 
planes. In ground-based astronomy this becomes a par- 
ticularly difficult task, due to the variable atmospheric 
conditions which make it difficult to detect the relatively 
low polarization degrees which characterize most astro- 
nomi cal sources (a few percent; see for example iLerovl 
2000). These fluctuations, in fact, introduce flux varia- 
tions among different polarization directions which can 
be eventually mistaken for genuine polarization effects. 

This problem has been solved in a number of different 
ways reviewed bv lTinbergenl l)1996[) and to which we refer 
the reader for a detailed description. In this paper we will 
focus on the so-called dual-beam configuration, which is 
the most popular one for instruments currently mounted 
at large telescopes. Despite new technologies, the basic 
concept of astronomical dual-beam polarimeters (see for 
example lApoenzellerl Il967t iScarrot et alJ Il98.1t) has re- 
mained unchanged. A mask placed on the focal plane, 
which prevents image (or spectra) overlap, is followed 
by a Wollaston prism, which splits the incoming beam 
into two rays characterized by orthogonal polarization 
states and separated by a suitable angular throw. The 
rotation of polarization plane is usually achieved with 
the introduction of a turnable retarder plate (half or 
quarter wave for linear and circular polarization, respec- 
tively) just before the Wollaston p rism (see for example 
iSchmidt. Stockman fc Smith! 1 1992ft . Recently new solu- 
tions have been proposed, in order to fully solve the prob- 
lem in one single exposure (see Oliva 1997 and Pernechele 
et al. 2003 for an example application), but so far they 
have been implemented in a few cases only. 

Alternatives to Wollaston-based systems have been de- 
vised. They are mainly based on the charge transfer 

*This paper is partially based on observations made with ESO 
Telescopes at Paranal Observatory under programme IDs 066. A- 
0397, 69.C-0579, 069.D-0461 and 072.A-0025. 
Electronic address: fpatat@cso.org, mromanie@eso.org 



in CCDs, which allows an on-chip storage of two dif- 
ferent polarization states, which are obtained rotating 
a polarization modulator. After t he pioneering work of 
iMcLean. Aspin fc Reitsenial l|1983j) . this technique, orig- 
inally proposed by P. Stockman, h as been succe ssfully 
applied in a number of instruments l)McLeanlll997[L 

In this work we address the most relevant problems 
which are connected to two-beam polarimetric observa- 
tions and data reduction. The paper is organized as 
follows. In Sec. we introduce the basic concepts of 
the problem and in Sec. we recall the analytical ex- 
pressions for the uncertainties of polarization degree and 
angle, which are then compared to Monte-Carlo simu- 
lations in Sec. fU In the same section we also recap the 
basics of polarization bias. Sec.EJdeals with the effects of 
background on the polarization measurements and Sec. El 
treats the flat-fielding issues. SecEJis devoted to the de- 
viations of the Wollaston prism from the ideal behaviour, 
the consequences of retarder plate defects are addressed 
in Sec. |H1 while the effects of post-analyzer optics are dis- 
cussed in Sec. El Sec. EHis dedicated to the instrumen- 
tal polarization and Sec. ^2 deals with the case of VLT- 
FORS1. Finally, in Sec. El we discuss and summarize 
our results. 

2. BASIC CONCEPTS 
The polarization state of the incoming light can be de- 
scribed t hrough a Stokes vecto r S(I, Q, U, V) (see, for 
example, iChandrasekharl Il9 50) . Its components, also 
known as Stokes parameters, have the following mean- 
ing: / is the intensity, Q and U describe the linear po- 
larization and V is the circular polarization. Linear po- 
larization degree P and polarization angle \ are related 
to the Stokes parameters as follows: 
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\/Q 2 + U 2 



Q 2 + U 2 
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where we have introduced the normalized Stokes pa- 
rameters Q = Q/I and U = U/I. The above relations 
can be easily inverted to yield: 



Q = P cos 2 X ; U = P sin 2 X 



(3) 



Finally, the circular polarization degree, not discussed 
in this paper, is simply P c — V = V/I. For the sake of 
clarity, we will set V—0 and neglect all circular polariza- 
tion effects throughout the paper. 

The ideal measurement system for linear polarization 
is composed of a half-wave retarder plate (HWP) fol- 
lowed by the analyzer, which is a Wollaston prism (WP) 
producing two beams with orthogonal polarization direc- 
tions. In general, each of these elements can be treated 
as a mathematical operator that acts on the input Stokes 
vector S (see for examDle lShurcliHl962tlGoldsteinl 2003'). 
What one measures on the detector is the intensities in 
the ordinary and extraordinary beams at a given HWP 
angle 9i, which are related to the Stokes parameters by: 



fo,i = W + Q cos46». i + C7 sin 40*] 
fs,i = k[I-Q cos40 i -E7 sin 40,] 



(4) 



If the observations are carried out using TV positions for 
the HWP, the whole problem of computing I, Q and U 
reduces to the solution of the 2N linear equations system 
given by Eqs.0] It is clear that, having three unknowns 
(I, Q and U), at least N=2 HWP position angles have 
to be used. 

Introducing the normalized flux differences Fi 



F, 



fo.i — Ie. 



fo.i + fE, 



(5) 



and noting that fo,i + Je,-. 
following N equations: 



/, Eqs. 0] reduce to the 



Fi = Q cos 40, + U sin 40< = P cos(40j - 2 X ) (6) 

We note that each F parameter is totally determined 
by a single observation, and it is therefore indepen- 
dent from sky conditions changes. It is also worth 
mentioning that alternative approaches to the normal- 
ized flux ratios exist. One examp le can be found in 
iMiller. Robinson fc Goodrich! l(l"987|h 

In principle, one can use any set of HWP angles to 
solve the problem, but is is easy to show that adopting 
a constant step A8=it/8 is the optimal choice. In fact, 
besides minimizing the errors of the Stokes parameters, 
this choice makes the solution of Eqs. trivial: 



2 7T 

Q = J7 Fi cos ( 2 

4=0 

2 JV_1 7T 
i=0 



(7) 



(8) 



Finall y, it prevents the "power leakage" (see, for in- 
stance, IPress et, alJ 119991) when one is to perform a 
Fourier analysis (sec below). 

In the ideal case, the normalized flux differences Fi 
obey Eq. which is a pure cosinusoid. Since all possible 



effects introduced by the HWP must reproduce after a 
full revolution, it is natural to consider them as harmon- 
ics of a fundamental function, whose period is 2tt. 

Therefore, if 0j — iri/8, Eq. EJcan be rewritten as the 
following Fourier series: 



Fi = ao 



N/2 2m 
J2 a k sin(fc— ; 

i=l 



h cos(fc^) 



where the Fourier coefficients are given by 



N—i 



ao 



a-k - 



i=0 



2 N ~ 1 2iri 
- ^ cos(fc— ) 



(9) 



i=0 



2 Ar_1 2m 

bh = — Fi sin(fc ) 

N ^ K N ' 

i=0 

which are valid for N = 4, 8, 12, 16. Comparing Eqs.|Hl 
with Eq. it is clear that the polarization signal is car- 
ried by the k = N/4 harmonic. In a quasi-ideal case, all 
Fourier coefficients are expected to be small compared 
with ajv/4 and 6^/4 and deviations from this behaviour 
could arise from a number of effects. For such an ap- 
proach to the error analysis and for the me aning of the 
variou s harmonics, the reader is referred to lFendt et alJ 
(1996). Here we just notice that the ao term, which 
should be rigorously null in the ideal case, is related to 
the deviations of the WP from the ideal behaviour (see 
Sec.EJ. 

In general, a Fourier analysis is meaningful when 
iV=16, and can reveal possible problems directly related 
to the HWP quality (cf. Sec. [HJ) ■ In most cases, though, 
due to practical reasons, one typically uses N— =4 and, in 
that circumstance, a different error treatment is required. 

3. ANALYTICAL ERROR ANALYSIS 

Under the assumption that all relevant quantities are 
distributed according to Gaussian laws, one can analyti- 
cally derive simple expressions for the corresponding er- 
rors of the final results. As we will see in the next sec- 
tion, this assumption is not always correct and, when 
this happens, a numerical treatment is required in order 
to test the analytical results and their range of validity. 
Assuming that the background level is the same in the 
ordinary and extraordinary beams and that the read out 
noise can be neglected, the analytical expression for the 
abs olute error of P can be readily deriv ed (see for exam- 
ple |MiiIi^^^binso^^^oodric3^83) propagating the 
various errors through the relevant equations 1 : 



crp 



^/N/2SNR 



(10) 



where SNR is the signal-to-noise ratio of the inten- 
sity image (fo + /s)- The signal-to-noise one expects 

1 Here we consider photon shot noise as the only source of ran- 
dom errors. Another potential source is represented by the mis- 
positioning of the HWP with respect to the optimal angles. How- 
ever, as analytical solutions and numerical simulations show, with 
the typical positioning accuracy nowadays attainable (< 1°), the 
associated error of the polarization degree and angle is negligible. 
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to achieve in the polarization degree, SNRp = P/ap, is 
simply given by: 

'N 



SNR P = \ — P SNR 
\i 2 



As for the error of x, this is given by: 



1 



(Tp 



2y/N/2 P SNR 2P 



(11) 



from which it is clear that, at variance with the po- 
larization degree, the accuracy of the polarization angle 
does depend on the intrinsic polarization degree. 

4. MONTE CARLO SIMULATIONS 

The analytical treatment presented in Sec. |3] relies on 
the assumption that all relevant variables obey Gaussian 
statistics. Numerical simulations are required in order to 
derive more realistic distributions and to verify the va- 
lidity of the analytical results. One can easily implement 
the concepts we have developed until now in a Monte 
Carlo (MC) code, which also allows higher sophistication, 
like the inclusion of Poissonian noise. With this tool one 
can readily investigate the effects of non-Gaussian dis- 
tributions of the derived quantities, the most important 
of which is the systemat ic error of t he po larization, as it 
was first pointed out bv iSerkowskil l)1958j) . 

4.1. Linear Polarization Bias 

Due to the various noise sources, the vector compo- 
nents Q and U are Normally distributed, but since P 
is defined as the quadrature sum of Q and U, the sta- 
tistical errors always add in the positive direction, lead- 
ing to a systematic increase of the estimated polariza- 
tion degree, thus introducing a bias. The problem was 
addressed by sever al authors, both with analytical and 
nume rical methods iSerkowskilQ SStlWardle fc Kronberd 
T97llSimmons fc Stewartlll985tlClarke fc Stewartlll986t 
" Sparks fc Axonlll999|h We refer the reader to those pa- 
pers for a detailed description of the problem, while here 
we just recall the basic concepts and we apply them to 
our case. 

The polarization bias is usually quantified using a ro- 
bust estimator, supposedly giving a statistically signifi- 
cant representation of the observed value, which is then 
compared with the input polarization in order to derive 
the systemat ic correction. Differe nt choices have been 
adopted (see ISparks fc Axon! 1 19991) Follo wing the con- 
siderations bv TWardfe^TlCronbergl l)1974|) we have here 
adopted the mode (P) of the distribution in order to esti- 
mate the bias, that we therefore define as AP — (P)—Pq, 
where Pq is the input polarization. Once applied to to 
the observed data, the bias correction AP tends to re- 
store the symmetry of the dey iation distribution (see for 
example iSparks fc Axonlfl999l their Fig. 4). 

In Fig. we show the results of our MC simulations 
for the estimated RMS error of the polarization o~p (up- 
per panel) and p olarization bias AP (l ower panel) for 
7V=4. Following ISparks fc A'xonl l)1999|) . we have used 
r\ = Pq ■ SNR and (rj) = (P) • SNR as independent vari- 
ables in our plots. As we had anticipated, ap follows 
the analytical prediction of Eg. 1 101 when r\ >2. For lower 
values of rj, ap tends to be systematically smaller than 




OA 0_6 0.S 
log (N/2)"> <T7> 

Fig. 1. — Upper panel: comparison between the RMS error on 
the polarization degree from MC simulations (circles) and Eq. I1UI 
The dashed line traces the expected RMS error for the Rayleigh 
distribution (see text). Lower panel: bias estimated using the mode 
(circles) a nd the average (crosses). Fo r comparison, the solid curve 
traces the Wardlc & Kronbcrg ( 1974) solution for the mode, while 
the dashed line is the TsTarks&liSonl 119991) solution for the aver- 
age. 



the analytical prediction and it converges to the value ex- 
pected for the Rayleigh distribution (dashed line), which 
becomes a very good approximation for r] <0.5, provided 
that SNR >3. For intermediate valu es of ij, the distribu- 
tion is described by a Rice function ( R.ice lll944j) . In con- 
clusion, one can safely use the analytical solution given 
by Eq.Unifor r\ >2 only. 

As for the polarization bias, we have plotted it as a 
function of measurable quantities, namely the signal-to- 
noise and the observed polarization (P). 

The results of our MC simulations, as plotted in the 
lower panel of Fig. ^ are i n good agreement w i th the 
analytical solution found bv lWardle fc Kronberd (|1974D 
for the same statistical estimator. For comparison, we 
have also plotted the results one obtains when the av- 
erage is adopted (crosses). For \/N/2 (77) >4 the re- 
lation between log(AP/(P)) and \og(y/W/2 (77)) is well 
approximated by a linear law. A least squares fit gives 
the following result: 



0.62 



1.92' 



which can be used to correct the observed polarization 
values according to the input signal-to-noise ratio, mea- 
sured polarization and number of HWP positions. In 
general, the bias effect is present even at reasonably high 
values of SNR when the polarization is small and ap 
and AP tend to become similar, so that the systematic 
bias correction is comparable to the random uncertainty 
of the polarization. This is better seen in Fig. [3 where 
we have plotted the ratio AP over up as a function of 
Po/o~p = SNRp deduced from our simulations. As an- 
ticipated, for low values of SNRp, the ratio between AP 
and ap tends to unity, with some variations among differ- 
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Fig. 2. — Comparison between systematic bias AP and random 
error op for different estimators as a function of input polarization 
signal-to-noise ratio. 

ent estimators. For SNRp >3 the bias correction is less 
than 10% of the expected accuracy, and it is, therefore, 
negligible. Moreover, above that threshold, all estima- 
tors give practically identical results. 

In Fig.Qlw e have plotted, for comparis on, the function 
computed bv lSimmons fc Steward 1)19851) . who have used 
a Maximum Likelihood estimator in order to evaluate 
the bias. As these authors have shown, this is the best 
estimator for SNR p <0.7, while for SNR P >0.7 the 
mode, first used bv iWa rdle & Kron berel l)1974[L should 
be used. 

5. THE EFFECTS OF THE BACKGROUND 

Until now we have assumed that one is able to perfectly 
subtract the background contribution. This is most likely 
the case when one is to perform polarimctric measure- 
ments on point-like sources, since in that situation local 
background subtraction is in most cases straightforward. 

We remark that the background, whatever its nature 
is, must be subtracted before t he calculation of normal- 
ized Stokes parameters (see alsolTinbcrgc nll996() , so that 
possible background polarization is vectorially removed. 

If we assume that the object is characterized by P 
and Xo and the background by Pb and Xb, the- two 
polarization fields can be expressed using the Stokes 
vectors defined as S (I ,I P cos 2x ,7 P sin 2% ) and 

Sb(h, hPb cos 2\b, hPb sin 2xb), where we have neglected 
any circular polarizat ion. Since Stokes vect ors are addi- 
tive (see for example lChand rasckharl ll950jL the result- 
ing polarization field is described by S = S + Sb and, 
therefore, the total polarization is given by the following 
formula: 

P = 7^7" Vl + r 2 + 2rcos[2( Xo - *„)] 

where r — (IbPb)/ \I Po), i-e- the ratio between the 
polarized fluxes of background and object. The corre- 
sponding polarization angle is 

1 sin 2xo + r sin 2xb 
X = — arctan 

2 cos 2xo + T cos2xb 



Clearly, the background is going to influence signifi- 
cantly the object when r > 1. For r ~ 1 one can write: 



P~ + cos[2( Xo - 

which implies that, for comparable polarized fluxes, 
the resulting polarization is nulled when the polarization 
fields are perpendicular (\xo ~ Xb\ = 7r /2)- 

6. FLAT FIELDING 

One of the basic problems one has to face when reduc- 
ing the data produced by dual-beam instruments is the 
flat-fielding. Due to the fact that image splitting occurs 
after the focal mask, the collimator and the HWP, one 
would in principle need to obtain flat exposures with all 
optical components in the light path. Unfortunately, in 
all practical conditions, this introduces strong artificial 
effects due to the strong polarization typical of flat-field 
sources (either twilight sky or internal screens). In prin- 
ciple one can reduce this effect using the continuous rota- 
tion of the HWP as a depolarizer. This is implemented, 
for example, in EFOSC2, cu rrently mounted at the ESO- 
3.6m telescope (Patat 1999) and it is effective only if the 
HWP rotation time is much shorter than the required ex- 
posure time. The depolarizing effect can also be achieved 
by averaging flats taken the same set of HWP angles used 
for the scientific exposures. In fact, with the usage of the 
optimal angle set (see Sec.|2J), one has that 



= 




and a similar expression for fs,u which do not contain 
any polarization information. The problem is that this 
is true only if the source is stable in intensity, which is 
surely not the case for the twilight sky and probably not 
really true for most lamps. 

An alternative solution (at least for imaging) is the us- 
age of a set of twilight flats obtained without HWP and 
WP. If on the one hand this eliminates source polariza- 
tion, on the other hand it does not allow for a proper 
flat field correction. In fact, while the pixel-to-pixel vari- 
ations are properly taken into account, the large scale 
patterns are not, due to the beam split that maps a given 
focal plane area into two distinct regions of the post- 
WP optics and detector. Moreover, these calibrations do 
not carry any information about possible spatial effects 
introduced by the HWP and the WP. However, as the 
simulations show, this problem becomes milder if some 
redundancy is introduced. For example, if one uses N=A 
HWP positions, the ordinary and extraordinary ray will 
just swap when the angle differs by 7r/4 within each of 
the two redundant pairs (/q = —f^ 2 , see Eqs.0J. This 
tends to cancel out the flat-field effect and becomes more 
efficient if the maximum redundancy (AT=16) is used. It 
must be noticed, however, that time dependent effects, 
like fringing, may affect the redundant pairs in a different 
way, therefore decreasing the cancellation efficiency. 

7. EFFECTS OF A NON-IDEAL WOLLASTON PRISM 

So far we have assumed that our system is described 
by Eqs. 01 i.e. that the Wollaston prism splits incoming 
unpolarized light into identical fractions. A deviation 
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from this ideal behaviour can be described by the intro- 
duction of a new parameter t in Eqs. 0] which can be 
reformulated as follows: 



SNR = 1607 x=922.4 y= 7 1 9. 7 r=1.1 arcmin 



fo,% = t[I + Q cas49i + U sin 40* ] 
fs,i = (l-t)[I~Q cos 46», - U sin 4^; 



(12) 



An ideal system is obtained for t = h. Now, for 
unpolarized light (Q—U=0), these new equations give 
fo,i = tl and fE,i = (1 — t)I for all HWP angles, so that 
all normalized flux differences turn out to be identical, 
i.e. Fi = 2t — 1. Therefore, the value of t can be directly 
estimated observing an unpolarized source. 

In the simplest situation, where N—2, neglecting the 
presence of the t term would lead to a spurious polariza- 
tion degree P — \/2 (2t — 1), with a polarization angle 
X = 7r/8. It is interesting to note that this is not the case, 
for example, when N=A. In that situation, in fact, be- 
cause all redundant F's are identical, Eqs. [7| and |S] would 
correctly yield null Stokes parameters. 

The problem when the incoming light is polarized is 
more complicated, since the normalized flux differences 
are not anymore a linear combination of Q and U, as one 
can verify from Eas. 1121 
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Fig. 3. — Example of Fourier Analysis applied to archival VLT- 
FORS1 observations of a bright star in the V passband (see text). 
Upper panel: normalized flux differences. Partial reconstructions 
using 8 harmonics (solid curve) and the 4-th harmonic only (dashed 
curve) are traced. The dashed horizontal line is placed at the 
average of F values (ao). Central panel: residuals after subtracting 
the fc=l,2,4 components. Lower panel: harmonics power spectrum. 



Fi = 



k + Qcos 49i + U sin 40j 
1 + kQ cos 40i + kU sin 46 l 



(13) 



where we have set k = It — 1 (|re| <1, re=0 in the ideal 
case). If re is known, one can correct the observed fo and 
Je dividing them by 2t and 2(1 — t) respectively before 
following the procedure adopted for the ideal case. Of 
course, if the source has a known polarization (e.g. a po- 
larized standard), one can use this information together 
with the observed F ratios to derive re for each HWP 
angle, according to the following relation: 

Q cos 40; + U sin 40 i - Fi 



F l (Qcos46»i + £/sin40 i 



If the input polarization is unknown, then one can in 
principle derive re from the observations themselves, pro- 
vided that N >4. In fact, after introducing the parame- 
ter 



9j 



1 

~F 



F 



3+1 



(14) 



and using Eqs. 1131 it is easy to demonstrate that 



= m ± \ 9 



1 



(15) 



where j = 1, ...,N/2 — 2 and the positive sign refers 
to the case gj < 0. For instance, for N = 4, one can 
determine two independent estimates of re, which can be 
averaged to improve on the accuracy. 

It is interesting to note that,_when P -C l, Eg. 1131 can 
be approximated as Fi ~ re + Q cos 46^ + U sm 48i. If N 
is a multiple of 4 (i.e. if the F function is sampled for an 
integer number of periods 7r/2), one has that 



N-l 

N ^ 

i=0 



ao 



(16) 



which clarifies the meaning of the ao term in the 
Fourier series (see Eq. |5J). 

It is worth mentioning that the redundancy on the F 
parameters does eliminate the effects of a non ideal WP 
to a large extent. For example, a blind application of 
Eqs.|7|and|21to the case of N=4 gives the following result: 



2 

N 



N-l 

F t cos(40i) 



1 



»=i 



1 - re 2 Q 2 



and a similar expression for U . Jf the polarization is 
small (Po < 0.1), we have that re 2 Q 2 -C 1 and therefore 
the application of the ideal case procedure to a non ideal 
situation would lead to a value of Q which is (1— re 2 ) times 
smaller than the real one. Since the same is true for U, 
the resulting polarization P will also be (1 — re 2 ) times 
smaller than the input value, while the polarization angle 
remains unchanged. For example, if |re| < 0.2, ap/P is 
less than 4%. 

Eas. 1121 describe a particular case only, where the in- 
coming unpolarized flux is distributed into two fractions 
t and (1 — t). More in general, one should replace the 
term (1 — t) with an independent parameter s, so that 
the fraction of light split by the WP in the ordinary and 
extraordinary ray become uncorrelated. Using the same 
procedure, it is easy to demonstrate that one can esti- 
mate the ratio Kj = (t — s) /(t + s) still using Eq. ^| 

Another effect we have investigated is the possibility 
that the difference in polarization direction between the 
ordinary and the extraordinary ray of the WP is not tt/2. 
If we call Aa the deviation from this ideal angle, using 
the general expression of t he Mueller matr ix for a linear 
polarizer (see for example Goldstein 2003() and deriving 
the expressions for the normalized flux ratios, one gets: 

Fi ~ ^[cos(46' l -2Aa)+cos46i l ] + Y[ sin ( 46, «- 2AQ! )+ sin46, l 
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where the approximation is valid under the assump- 
tion that ? « 1. With the aid of this expression it is 
easy to conclude that for reasonably small values of Aa 
(<10°), the implied errors are of the order of 0.05% on 
the polarization degree and 5° on the polarization angle, 
irrespective of the number of HWP positions used. 

8. HWP DEFECTS 

In the ideal case, the normalized flux differences are 
modulated by the HWP rotation according to Eq. 
which is a pure cosinusoid. If defects are present on the 
HWP, like dirt or inhomogeneously distributed dust, one 
can expect spurious flux modulations, which are not re- 
lated to the polarization of the incoming light and can 
reduce the performance of the instrument. As a conse- 
quence, error estimates based on pure photon statistics 
are systematically smaller than the actual ones. 

These kinds of problems can be investigated with the 
aid of Fourier Analysis, following the procedure we have 
outlined at the end of Sec. El This method becomes par- 
ticularly effective when the observations are taken sam- 
pling the full HWP angle range, i.e. 2n which, given the 
choice of the optimal angle set 9i=iri/8, implies A=16 
retarder plate positions. In these circumstances, one is 
able to determine the Fourier coefficients and bk for 
the first 8 harmonics, the fundamental (k = 1) being 
related to local transparency fluctuations which repeat 
themselves after a full revolution, like dirt or dust. By 
definition (cf. Eq. |SJ, the fourth harmonicjs directly re- 
lated to the linear polarization, i.e. a^ = Q and 64 = U. 
All other harmonics, with the only remarkable exception 
of the second one (k=2), are simply overtones of harmon- 
ics with lower frequencies and include part of the noise 
generated by the photon statistics, which is present at all 
frequencies and is therefore indicated as white noise. For 
this reason, the global random error is often estimated as 
the signal c arried by the harm onics having fc=3,5-8 (see 
for example iFendt et alJll99fl their Appendix A. 3). 

The second harmonic deserves a separate discussion. 
Ideally, the HWP operates as a pure rotator of the input 
Stokes vector, with the advantage that one does not need 
to rotate the whole instrument in order to analyze differ- 
ent polarization planes. In the real case, being the HWP 
usually constructed using bi-refringent materials, it is af- 
fected by the so-called pleochroism. This is a wavelength 
dependent variation of the transmission that takes place 
when the direction of the incoming light is changed with 
respect to the crystal lattice. Due to the way the HWP 
is manufactured, the crystals have an axial symmetry, 
which gives a period of n. Therefore, this effect is seen 
as the k=2 component. 

In Fig. |31 we show a real case, where we have applied 
this analysis to archival data obtained with the FOcal Re- 
ducer/low dispersion Spectrograph (hereafter FORS1), 
which is currently mount ed at the Cas segrain focus of 
ESO-VLT 8.2m telescope (|Szeiferdl200^> . 

A bright (and supposedly unpolarized) star was ob- 
served using 7V=16 HWP positions. First of all, the 
Fourier Analysis indicates the presence of a small de- 
viation of the WP from the ideal behaviour (k « 
cto ~4.1xl0 -3 , see Sec. EJ). Then, a clear polarization 
is detected, at the level of about 0.4% while all other 
components are smaller than 0.05% (this polarization is 
actually an instrumental effect present in FORS1. See 



Sec. HOpi . The effective significance of harmonics other 
than fc=4 can be judged on the basis of the expected er- 
rors of the Fourier coefficients. For example, using the 
expression of a>k one finds that 



2 A , , 2m 

<7„, = — . > COS z K 

a " N SNR N 

With this kind of analysis, one can see that in the 
example of Fig. |31 and are consistent with a null 
value for fc=3,5-7 (see the central panel). As for the 
fc=l,2 harmonics, the Fourier coefficients are non null at 
a 2a level. Since for the test star it was SNR ^1600, it 
is clear that to detect k—1,2 harmonics of this amplitude 
(0.05%), a SNR >3000 is required. 

It is important to notice that the presence of these non- 
null harmonics is implicitly corrected for when one has 
a sufficient number of HWP positions covering the max- 
imum period 2tt. In the most common case, where N=A 
angles spaced by tt/8 are used, one can derive the funda- 
mental (i.e. linear polarization, period n/2) and the first 
overtone (period 7r/4) only. The latter corresponds to 
the fc=8 component of the N—16 cases, which therefore 
carries the high frequency information only. As a con- 
sequence, if other harmonics are present, they are not 
properly removed and contribute to the final error, prac- 
tically setting the maximum accuracy one can achieve, ir- 
respective of the SNR. Numerical simulations performed 
assuming a virtually infinite SNR show that in the pres- 
ence of k = 1 and k — 2 components, the usage of 7V=4 
HWP angles leads to systematic errors which are of the 
same order of the amplitude of the two harmonics. From 
this and the example reported in Fig. |21 one can es- 
timate that the absolute maximum accuracy reachable 
with FORS1 using N=4 is of the order of 0.05%. 

Another typical problem which affects the retarder 
plates is the chromatic dependence of the angle zero 
point. This is usually measured by means of a Glan- 
Thompson prism and it can change by more than 5° 
across the optical wavelength range. The computed po- 
larization angle can be corrected simply adding the HWP 
angle offset for the relevant wavelength (or effective wave- 
length in the c a se of broad band imaging). See, for ex- 
ample, Hiiferll (Hop). 

Finally, we have investigated the effects produced by a 
deviation A/3 from the nominal phase retardance of an 
HWP (it). Using the general expressi on of the Muelle r 
matrix for a retarder (see for example IGoldsteml [2003) . 
the normalized flux ratios turn out to be 

F i = Q [cos 40* - sin 2 20^1- cos A/3)] + 

+ ~Usin40i (1 + cos A/3) + V sin 20* sin A/3 

from which it is clear that the measured linear polar- 
ization depends on the circular polarization of the input 
signal. For V=0 and A/3 <10° the corresponding abso- 
lute error of the polarization degree is less than 0.05%, 
while the outcome on the polarization angle is negligible. 
For V 7^0 the exact effect depends on the ratio between 
circular and linear polarization degrees. For example, for 
A =10° and Q = U = V =0.01, the absolute error of the 
computed polarization degree is about 0.1% for iV=4, 
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which decreases to 0.01% for N—16. It is worthwhile 
noting that this defect would be detected by a Fourier 
analysis as a component with a period 7r and whose in- 
tensity is |V"sin A|. 

9. EFFECTS OF POST-ANALYZER OPTICS 

Typically the analyzer is followed by additional optics, 
like filters, grisms and camera lenses which, due to their 
possible tilt with respect to the optical axis, may be- 
have as poor linear polarizers. In the most probable case 
where the polarization is produced by transmission (see 
also next section), the properties of post-analyzer (PA) 
components can be described by the following approxi- 
mate Mueller matrix: 

(I BcBs 0\ 
MpA ~ A \B S CO 

V o o o cj 

where A s»l, C sal, c = cos2</?, s = sin2<p and (p is 
the polarization angle (which can change across the field 
of view) while B is related to the polarization degree 
introduced by the PA optics. This expre ssion can be 
deduced from the general formulation (see lKellerll2002l 
Eq. 4.63) after applying the usual matrix rotation (see 
for example lKeTIerl2002l Eq. 2.5). If S = (I , Q , U Q , V ) 
is the input Stokes vector, the effect of PA optics can be 
evaluated computing the Stokes vectors that correspond 
to the ordinary and extraordinary beams produced by 
the WP, transforming them using the operator Mpa and 
using the resulting intensity components to compute the 
normalized flux differences Fi. After simple calculations, 
one arrives at the following expression: 

Fi ~ B cos 2ip + Q cos 40, + Uo sin 0j 

were we have assumed that \B\ <C 1, i.e. that the 
linear polarization induced by the PA optics is small. 
Given this expression, it is clear that the redundancy in 
the HWP positions (N=4, 8, 16) eliminates this prob- 
lem, since the additive term B cos 2ip is not modulated 
by the HWP rotation, while for A^=2 the derived polar- 
ization degree and angle would be affected, possibly in 
a severe way. If the optimal HWP angle set has been 
used, it is easy to verify that Bcos2ip = X^o* Fi/N, 
which is identical to Eq. El This means that, in a first 
approximation, it is not possible to disentangle between 
an imperfect WP and the presence of polarization in 
the PA optics. Therefore, the fact that ao ~0.4% in 
Fig. |21 can actually be attributed to both kinds of prob- 
lems. The PA optics effect becomes definitely stronger 
when these include highly tilted components, like the 
grisms. This is very well illustrated by the two examples 
of Fig. 01 where we show the results obtained using VLT- 
FORS1 archival data of a highly polarized star (Vela 1 
95, a=09:06:00, £=-47:19:00) and an unpolarized star 
(WD 1615-154, a=16:17:55, <5=-15:35:51) 2 , which were 
observed on the optical axis, where the instr umental po- 
larization is known to be null l)Szeifertll20 02) . 

In both cases the polarization degree deduced using 
N—2 (central panels) is markedly different from that 

2 See http://www.eso.org/instruments/fors/inst/pola.html 
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Fig. 4.— VLT-FORS1 observations of Vela 1 95 (left panels) 
and WD 1615-154 (right panels). The plots show the linear po- 
larization derived with N=A HWP positions (upper panels), N=2 

(middle panels) and ^2fLg Fi/N f° r N=A (lower panels. See text 
for more details). The original spectra were obtained with the 300V 
grism and a slit of 1"; for presentation they have been binned to 
25 A. The empty circles in the upper left panel mark the broad- 
band polarimetric measurements for Vela 1 95 (UBVRI, from left 
to right). 



derived with A^=4 (upper panels) and the deviation is 
particularly severe for the unpolarized object. As one 
can finally notice, the resulting values of B cos 2(p show a 
strong wavelength dependency and are higher than 5% at 
about 800 nm. It is interesting to note that B cos 2ip ^0 
at about 450 nm, i.e. at the wavelength where the anti- 
reflection coatings are optimized (see next section) . This 
fact, together with the marked wavelength dependency 
and the much lower level seen in broad band imaging (see 
Fig. strongly suggest that the effect seen in Fig. 01 is 
indeed produced by the tilted surfaces of the grism. 

10. EFFECTS OF INSTRUMENTAL POLARIZATION 

So far, we have assumed that all optics preceding the 
analyzer do not introduce any polarization. Of course, 
this is not generally true (see for example iTinbergenl 
119961 and Leroy 2000 for a general introduction to the 
subject). 

To show the effect of instrumental polarization, we as- 
sume that the pre-analyzer optics, which include tele- 
scope mirrors, collimator, HWP and so on, introduce an 
artificial polarization, which depends on the position in 
the field. For the sake of simplicity, we assume these op- 
tics act as a non-perfect linear polarizer, characterized 
by a position dependent polarization degree p(x, y) and 
polarization angle tf(x, y). This can be described by the 
following Mueller matrix: 

1 I 1 pc ps \ 

Mi(x, y) = I pc 1 — ps 2 psc 

\ps pes 1— pc 2 / 

where we have neglected circular polarization and we 
have set s = sin 2tp and c = cos 2<p. For p=0 one obtains a 
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totally transparent component, while p=l gives an ideal 
linear polarizer. 

If ^(/o, Qo, Uq) is the Stokes vector describing the 
input polarization state, it will be transformed by pre- 
analyzer optics into the vector Si = Mi ■ Sq before en- 
tering the analyzer: 

h (1 + p) = I + pcQo + psU 
Qi(l+p)=Qo+pcI Q -ps 2 Qo+pcsUo (17) 
Ui (1 + p) = U Q + psI Q + pcsQ Q - pc 2 U Q 

and, therefore, the measurements would lead to Si 
which would then need to be corrected for the instru- 
mental effect inverting Eas. 1171 provided that p(x, y) and 
ip(x, y) are known. 

Of course, if the observed source is known to be unpo- 
larized, p and ip can be derived immediately, for example 
by placing a single target on different positions of the 
field of view or observing an unpolarized stellar field. 

If the source is polarized, the problem becomes much 
more complicated, since Eqs. 1171 are non linear in c and 
s. The solution can be simplified assuming that p« 1 
and Pq <C 1, which is a reasonable hypothesis in most 
real cases, since instrumental polarization is typically less 
than a few percent. In these circumstances, Eas.1171 can 
be rewritten as follows: 

1 1 ~J„ 

Qi~Qo+pIo cos 2ip (18) 
Ui ~ Uo + plo sin 2ip 

It is important to note that the instrumental polariza- 
tion is not removed by the local background subtraction. 
Moreover, it is independent of the object's intensity; in 
fact, using the previous expressions one can verify that 

P=^Pi+P 2 + 2PoP cos[2( Xo - tp)] 

where Po and xo & r e the input polarization degree and 
angle. From this expression it is clear that when Po 2* 
p it is also P w Po, while in the case that object and 
instrumental polarization are comparable (p ss Po), the 
observed polarization is approximately given by 

P =s >/2fVl + cos[2(xo-¥0] 

which, according to the value of (xo — V 3 ), gives val- 
ues that range from to 2Po. It is important to notice 
that the main difference between instrumental polariza- 
tion and a polarized background is that the latter is ef- 
fective only when BGR > I (see Sec.|S}, while the former 
acts regardless of the object intensity and what counts is 
its polarization. 

With the aid of these approximate expressions 
fEas. IT8*|) . one can easily evaluate the instrumental po- 
larization, provided that the input polarization field is 
known and the observed source covers a large fraction of 
the instrument field of view. In fact, solving Eas. 1181 for 
if and p yields: 



and 

Qi -Qo Ui- Uo 

Pi - 7 7T- 1 Pi - -j 

1 cos 2ip 1 sin 2ip 
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Fig. 5. — FORS1 instrumental polarization map in the B band. 
The contours trace 0.3%, 0.6% and 0.9% polarization levels. Co- 
ordinates, expressed in arcminutes, refer to the geometrical center 
of the detector. 
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Fig. 6. — Upper panel: FORS1 instrumental polarization ra- 
dial profile for the B band. Each point is the result of a 30x30 
px binning in the original images. Radius, expressed in arcmin- 
utes, is computed from the geometrical center of the detector. The 
thick line traces a linear least squares fit, while the thin line is the 
polarimetric ray- tracing prediction. Lower panel: instrumental po- 
larization angle as a function of pixel polar angle. The solid line is 
not a fit to the data, but rather has unit slope and zero intercept. 



where pi and P2 are two independent estimates of p, 
that can be averaged to increase the accuracy. 

As it is well known, the night sky shows a polarization 
which varie s according to the e cliptic and galactic coor- 
dinates fsee lLeinert et al1l!998l for an extensive review). 
It is mostly dominated by the zodiacal light polarization, 
which reaches its minimum, below a few p ercent, at the 
anti solar position ( Roa ch &: Gordonlll973j) . Since this is 
not expected to vary on the scales of a few arcminutes, in 
principle, relatively empty fields represent suitable tar- 
gets for panoramic polarization tests, provided that the 
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TABLE 1 
Data for the empty fields 

Filter a 8 A j3 A - A P sky Xsky 

hh:mm:ss dd:mm:ss ° % ° 

B 03:32:17 -27:44:24 -41.1 -45.1 128.1 6.45 +68.4 

V 13:58:03 -31:22:21 218.6 -18.1 -159.7 1.74 -53.8 

I 20:36:08 -13:06:39 308.0 +5.3 -110.0 1.34 +68.7 



signal-to-noise ratio per spatial resolution element is of 
the order of several thousands. 

11. THE CASE OF FORS1 AT ESO-VLT 

To show an example application, we have performed 
a test using real data obtained with FORS1 at the 
ESO-VLT. In this instrument, the polarimetric mode is 
achieved inserting in the beam a super-achromat ic HWP 
and a WP, which has a throw of about 22" (jSzeifertJ 
2002). 

We have identified in the ESO archive three sets of 
data obtained in rather empty fields in B, V and I pass- 
bands. Equatorial coordinates (a, 5), ecliptic longitude 
and latitude (A, /?), helio-ecliptic longitude (A — A ), sky 
polarization degree and angle (P s ky, Xsky) are reported 
in Table for the different fields. In all three cases the 
signal-to-noise ratio achieved on the sky background I s k y 
in the combined images is larger than SNR ~200 per 
pixel. With such a signal and for a typical 1% polariza- 
tion, the bias effect is expected to be small (see Fig. ^| 
and the RMS error of the polarization degree, according 
to Ea. 1101 is of the order of 0.3%, while the uncertainty 
of x is about 9° (see Eq. In order to further in- 

crease the accuracy and to allow for outlier rejection, we 
have computed a clipped average in 30 x 30 px bins which, 
given the FORSI detector scale (0".2 px^ 1 ), translates 
into an angular resolution of 6". 

Since the instrumental polarization on the optical axis, 
measured with unpolariz ed standard stars, is smaller 
than 0.03% <|Szeifertil20 021 . one can be confident that the 
sky background polarization field (P s ky,Xsky) measured 
close to that area is not affected by spurious effects (the 
values are reported in the last two columns of Table ^) . 
Therefore, we can easily compute p(x, y) and (p(x, y) 
using the method previously outlined, where Io=I s ky, 

Qo = I sky Psky COS 2Xsky and £/n — I s ky Psky SHI ^Xsky- 

The results of these calculations are presented in Figs.Jol 

and |5] With the remarkable exception of the B band, 
the instrumental polarization of FORSI shows a quasi- 
symmetric radial pattern. For example, for the V filter, 
the instrumental polarization remains below 0.1% within 

1 arcmin from the geometrical center of the detector, 
while it grows up to ~0.6% at 3 arcmin, to reach the 
maximum, i.e. ~1.4%, at the corners of the field of view. 

This is illustrated in the upper panels of Figs. H3 El and 
1101 where we have plotted the estimated instrumental 
polarization for each 30 x 30 px bin as a function of 
its average distance r from the center. The deviation 
from a perfect central symmetry is distinctly shown by 
the dispersion of the points, which is larger than the 
measurement error. Particularly marked is the case of 
the B band, which shows a strong azimuthal dependence 
and deserves a separate discussion (see next section). For 
the V band there is a systematic deviation from central 
symmetry for a polar angle a = arctan(y/x) between 
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Fig. 7. — Same as Fig. l5l for the V band. The dark segment 
marked by a circle in the lower right corner indicates the values 
obtained from an unpolarized standard star. 



10° and 80°. This region is probably disturbed by the 
presence of a saturated star and a reflection caused by 
the HWP, visible in the input images. Excluding these 
points, a linear least squares fit to the V data gives 

p(r) = 0.012r + 0.046r 2 + 0.002r 3 

where r is expressed in arcminutes (Fig. [HJ solid curve). 
The / band shows the smoothest behaviour and the ob- 
servations are described very well by the following poly- 
nomial: 

p(r) = -0.017r + 0.105r 2 - 0.006r 3 

The absolute RMS deviations shown by the data from 
the best fits are of the order of 0.05%, so that in these 
two passbands the spurious polarization can be corrected 
with an accuracy which is comparable to that dictated 
by the photon statistics. In both cases, but especially in 
the / band, the pattern is remarkably radial, as shown in 
the lower panels of Figs. 151 and llOl where we have plotted 
if as a function of polar angle a. 

In order to verify these results for the V filter, we 
have carried out a test observing an unpolarized stan- 
dard star placed in the lower right corner of the de- 
tector. Measured polarization was P=0.92±0.04% and 
X=— 48°±1°.4 while, according to the previous analysis, 
the expected instrumental polarization in that position 
is p=0.96% and cp = —51°. 9, which are in very good 
agreement with each other (see also Fig. lower right 
corner). 

Once the instrumental polarization is mapped, one can 
correct for it using the approximate Eqs. 1181 which hold 
when p and Po are small and only if the instrumental 
polarization is produced by a linear polarizer preceding 
the analyzer. 

11.1. The cause of instrumental polarization in FORSI 

As we have seen, the spurious polarization detected in 
FORSI in V and / passbands h as a clear centra l symme- 
try, it is null on the optical axis lSzeifertJ (|2002T> . shows a 
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Fig. 8. — Same as Fig. IHlfor the V band. Fig. 10. — Same as Fig.^Jfor the I band. 
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Fig. 9. — Same as Fig.|2]for the I band. 



radial pattern and grows with the distance from the opti- 
cal axis. All these facts suggest that this must be gener- 
ated by the optics which precede the analyzer, i.e. within 
the collimator. In fact, when a light beam enters an op- 
tical interface along a non-normal direction, the compo- 
nent of the transmitted beam perpendicular to the plane 
of incidence is at tenuated, according to Fresnel equations 
(see for example iBorn fc Wolf3ll9"8 0') . As a consequence, 
the emerging beam is linearly polarized in a direction 
which is parallel to the plane of incidence. If the surface 
is curved, as it is the case of lenses, the incidence an- 
gle quickly increases moving away from the optical axis 
and this, in turn, produces an increase in the induced 
polarization. The effect becomes more pronounced if the 
lens is strongly curved, i.e. if the curvature radius is 
comparable to its diameter. Of course, on the optical 
axis, the incidence angle is null, so that no polarization 
is produced. Therefore, at least from a qualitative point 



of view, the polarization induced by transmission has all 
the required features to explain the observed pattern. 

The polarization induced by transmission can be easily 
evaluated using the appropri ate expressio n for the corre- 
sponding Mueller matrix fsee lKellerl2 002. Eg. 4.63). For 
a typical refraction index n=1.5 and an incidence angle 
of 30°, refraction through an uncoated glass would pro- 
duce a polarization of about B—1.7 % per optical surface. 
This polarization is usually reduced in a drastic way (i.e. 
down to 0.1-0.2%) by anti-reflection (AR) coatings. Nev- 
ertheless, since the effect of multiple surfaces is roughly 
additive, in the presence of numerous and pretty curved 
lenses, residual polarization can be non negligible. An- 
other important aspect is that this mechanism has no ef- 
fect on circular polarization, in agreement with the fact 
that no instrume ntal circular polariza tion has been mea- 
sured in FORS1 (Bagnulo et al]l2000jK 

We have run polarization ray-tracing simulations, in- 
cluding telescope mirrors, collimator lenses and AR coat- 
ings. This kind of calculation allows one to describe in 
detail the optical system, taking into account partial po- 
larization cancellation produced by symmetries within 
the optical beams and the depolarizing effect of AR coat- 
ings. 

The standard resolution collimator of FORS1 contains 
3 lenses and a doublet, all treated with a single-laye r 
MgF2 quarter- wave AR coati ng at 450 nm l)Seifertll 994'l. 
The ray-tracing calculations l)Avilall2005l) show that in- 
deed the polarization induced by transmission is not to- 
tally removed by the AR coatings. For V and / filters, a 
best fit to the simulated data gives a radial dependence 
which is very similar to the results we have derived from 
the experimental data. The deviation is maximum at the 
edges of the field of view, where the ray-tracing model 
gives a polarization which is ~0.08% and ~0.05% smaller 
than what is actually observed in V and / respectively 
(see also Figs. |H1 and upper panel, thin curves). Pos- 
sible explanations for this small discrepancy are to be 
identified with imperfections in the AR coatings and with 
the effects of non-orthogonal incidence on the HWP. 
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In principle, since single-layer AR coatings are opti- 
mized for one specific wavelength (450 nm in the case 
of FORS1), the residual polarization is expected to be 
higher at other wavelengths. Simulations have been run 
in order to sample the wavelength range 400-900 nm and 
they show that the expected wavelength dependency can 
be very well approximated by the following linear rela- 
tion: 

-^L 0.02+ 1.73 x 10~ 3 A 
p(550) 

where A is expressed in nm. This relation predicts 
pretty accurately the ~40% relative increase we indeed 
see passing from V to the / passband and can, thus, be 
safely used to predict the effect in the R band. 

According to the simulations, one would expect that 
the spurious polarization in B is about 25% smaller than 
in V . But, as we have already mentioned, this passband 
shows a rather weird behaviour and does not conform 
to the model predictions. In fact, the polarization pat- 
tern strongly deviates from central symmetry, display- 
ing a marked azimuthal dependence (Fig. left panel). 
This becomes more evident looking at the radial profile 
presented in Fig. HO (upper panel): the purely radial de- 
pendence is clearly disturbed by an asymmetric field. In 
some directions the polarization field grows much faster 
than in others, producing a great spread in the observed 
data and, in most of the cases, the observed polariza- 
tion is larger than what is predicted by the polarization 
ray-tracing (thin solid curve) . The deviations from a cen- 
trally symmetric pattern reach up to 0.5% (see Fig. Ill|) . 
making the correction in the B band quite difficult, and 
certainly not feasible using simple smooth functions as in 
the cases of V and I. Rather, a much more accurate cor- 
rection can be obtained by interpolating the map of Fig-El 
at the required field position. We must remark, however, 
that a rigorous correction for this secondary effect will 
be possible only once its physical reason is identified and 
its mathematical description is formulated. 

We have tried to reproduce the observed behaviour in- 
troducing defects in the system, like a weak linear po- 
larization from the HWP and the presence of linear po- 
larization in the post analyzer optics. In both cases the 
effect is completely different from what we see in the B 
band. Therefore, the physical reason of this phenomenon 
is still unclear (see also the discussion in the next sec- 
tion). What we can say here is that the deviation from 
central symmetry is present, though to a much smaller 
extent, also in the V band. This is shown by the contours 
at constant polarization, which are clearly box-shaped 
(see Fig. 0, while in the I band they are practically 
circular (see Fig. EJ . The conclusion is that this addi- 
tional effect, whatever its origin is, becomes more severe 
at shorter wavelengths. 

We must notice that our method heavily relies on the 
assumptions that the instrumental polarization is null on 
the optical axis and the sky background polarization is 
constant across the field of view. As a matter of fact, 
the night sky polarization is not very well studied. The 
only extensiv e analysis we could find in the lite rature 
is the one by Wolstencroft & Bandcrmann (1973), who 
concluded that the polarization structure varies in scale 
from a few degrees to about 30°, i.e. on scales which are 
much larger than the field of view of FORS1 (6.'8 x 6.' 
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Fig. 11. — Residual field obtained subtracting the ray-tracing 
model from the observed polarization in the B band. 



8). In order to explain the deviations we observe from a 
centrally symmetric pattern, one would need a variation 
in the sky polarization of the order 0.5% on a scale of a 
few arcmin. Even though this seems to be quite a large 
gradient, in principle we cannot exclude it. Only further 
tests will clarify the nature of the effect we see in the B 
band. 

12. DISCUSSION AND CONCLUSIONS 

Dual-beam polarizers coupled to 2D arrays provide 
a tool to perform panoramic imaging polarimetry and 
multi-object spectro-polarimetry. In these instruments, 
the atmospheric fluctuations problem is solved by obtain- 
ing simultaneous measurements of two orthogonal polar- 
ization states. Of course, the use of a WP has also some 
drawbacks, as the flat-fielding issue discussed in Sec. 
With the only exception of this feature, data reduction 
and analysis are totally similar to other polarimetric sys- 
tems, as we have shown with both analytical and numer- 
ical approaches (Sees. |3 and 0J. 

When the targets to be studied are extended and cover 
a large fraction of the field-of-view, accurate background 
subtraction becomes an issue, whose effects we have in- 
vestigated in Sec. This is particularly important when 
the background is not the simple sky background, but 
it has a complicated structure. This is the case, for in- 
stance, for a faint supernova projected onto a galactic 
spiral arm. 

Another problem that may reduce the performance of 
a dual-beam polarimeter is the imperfect behaviour of 
the WP. In Sec. H we have discussed this issue and pre- 
sented a test to determine possible deviations from the 
ideal case. As an example, we have applied it to the 
FORS1 archive data we have described in Sec. EH Using 
an object-free region roughly in the center of the field- 
of-view we have used Eqs. El an d El to compute t (see 
Eq.HSJl, which turns out to be i=0.502±0.001, i.e. per- 
fectly compatible with the value derived from the Fourier 
analysis (Sec.|SJ|. As we have shown, the redundancy in- 
troduced by having N > 4 strongly reduces this problem, 
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even in the cases where t differs by about 10% from the 
ideal case (i=0.5). This is the case also for the pres- 
ence of linear polarization in the post-analyzer optics, 
whose effects are practically eliminated by the redun- 
dancy (Sec. EJ) • 

Finally, we have addressed the instrumental polariza- 
tion issue, described its consequences and proposed an 
easy test to detect any spurious effect with rather high 
accuracy f Sec. 110^ . As an example, we have applied it to 
archival FORS1 data and we have detected an instrumen- 
tal polarization pattern, roughly centrally symmetric (for 
V and /) and with a radial dependency. The presence of 
this spurious polarization affects all objects placed at dis- 
tances larger than 1.'5 from the optical axis with intrinsic 
polarizations of a few percent or less. The problem be- 
comes particularly severe when p ~ P: in that case the 
measured Stokes parameters can be severely wrong. For 
objects filling most of the field of view, there always will 
be regions affected by this problem. Moreover, the cor- 
rect sky background estimate, which is absolutely neces- 
sary to recover the intrinsic object field in the outer parts 
of the galaxy, becomes impossible, if the instrumental 
polarization is not taken care of properly. The spurious 
field must be removed before one is able to estimate the 
background contribution. Both our data and ray-tracing 
simulations show that the effect is wavelength dependent. 
In the case of FORS1, a strong deviation from central 
symmetry is seen in the B band and we have interpreted 
this as a signature of an additional effect, yet to be ex- 
plained, not included in the ray-tracing simulations that, 
in contrast, reproduce quite well the observed data in V 
and /. 

One possible source of asymmetric instrumental polar- 
ization is the unrelieved stress birefringence in the opti- 
cal glasses, due to t hermal strain and mechanica l load- 
ing (see for example iTheocaris fc Gdoutosl ll979'). This 
phenomenon is known to introduce a retardance which, 
in turn, can change the polarization status of incoming 
polarized light (the effect is null if the light is unpolar- 
ized). Since the incoming radiation is certainly polarized 



Apponzeller, I., 1967, PASP, 79, 136 

Avila, G. 2005, ESO Technical Note INS-05/01 

Bagnulo, S., Szeifert, T., Wade, G. A., Landstreet, J.D. & Mathys, 

G. 2000, A&A, 389, 191 
Born, M. & Wolff, E. 1980, Principles of Optics, 6th ed., (New 

York: Pergamon Press) 
Chandrasekhar, S., 1950, Radiative Transfer, (Oxford: Oxford 

University Press) 
Clarke, D. & Stewart. B.G. 1986, Vistas Astron., 29, 27 
Fendt, C, Beck, R., Lesch, H. & Neninger, N. 1996, A&A, 308, 713 
Goldstein, D. 2003, Polarized Light, second edition, (New York: 

Marcel Dekker Inc.) 
Keller, C.U. 2002, in Astrophysical Spectropolarimetry, 

(Cambridge: Cambridge University Press, p. 303) 
Kendall, M. & Stuart, K. 1977, The Advanced Theory of Statistics, 

Vol. 1, (London: Charles Griffin & Co.) 
Leinert, Ch., Bowyer, S., Haikala, L.K., et al. 1998, A&AS, 127, 1 
Leroy, J.L. 2000, Polarization of Light and Astronomical 

Observation, (Singapore: Gordon & Breach Science Publishers) 
McLean, I.S., Aspin, C. & Reitsema, H. 1983, Nature, 304, 243M 
McLean, I.S. 1997, Electronic Imaging in Astronomy, John Wiley 

& Sons 

Miller, J.S., Robinson, L.B. & Goodrich, R.W. 1987, in 
Instrumentation for Ground-Based Optical Astronomy, Ed. L.B. 
Robinson, p. 157 

Oliva, E. 1997, A&AS, 13, 589 



by FORS1 in a differential way across the field of view, 
this would also imply that the secondary effect should 
be weaker where the centrally symmetric component is 
smaller. The fact that this is indeed the case (see Fig. lllfl 
and also that the retardance is expected to grow faster 
than A -1 seem to suggest that this is a plausible expla- 
nation for the asymmetric component. If this is indeed 
the case, then it is not possible to correct the measured 
linear polarization just subtracting vectorially the resid- 
ual field (like the one shown in Fig. Illfl simply because 
the effect of retardance depends on the polarization state 
of the incoming light. This requires a more sophisticated 
treatment, necessarily based on the exact knowledge of 
the physical mechanism and its mathematical description 
through Mueller matrices formalism. 

In general, instrumental polarization induced by trans- 
mission is most likely common to all focal reducers 
equipped with a polarimetric mode. While the overall 
pattern should be a general feature of these instruments, 
the exact radial dependence may change according to 
the optical design and the curvature of the lenses. The 
method we have described in this paper allows an accu- 
rate way of characterizing the instrument and a tool to 
correct for this effect. 

This paper is partially based on observations made 
with ESO Telescopes at Paranal Observatory under pro- 
gramme IDs 066.A-0397, 69.C-0579, 069.D-0461 and 
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